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Abstract 

Classical optical lithography is diffraction limited to writing features of a 
size A/2 or greater, where A is the optical wavelength. Using nonclassical photon- 
number states, entangled AT at a time, we show that it is possible to write features 
of minimum size AJ(2N) in an iV-photon absorbing substrate. This result allows 

one to write a factor of N 2 more elements on a semiconductor chip. A factor of 
N = 2 can be achieved easily with entangled photon pairs generated from optical 
parametric downconversion. It is shown how to write arbitrary 2D patterns by 
using this method. 
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Optical lithography has been the primary tool of the semiconductor 
industry for transferring circuit images onto substrates. However, diffraction 
effects in the typical masking approach limits the minimal resolvable feature size 

to the Rayleigh diffraction limit of A/2, where A is the optical wavelength. Hence, 
it has become necessary to use light of ever-shorter wavelengths to fabricate 
smaller features. Current production technology writes 180-220 nm features 
using KrF excimer laser light at 248 nm. New technological approaches consider 
light in the vacuum ultraviolet or soft x ray regime in order to obtain features at 
100 nm or below [1]. In all cases, the light is treated classically or, equivalently, as 
a stream of uncorrelated photons — an approach that leads to the Rayleigh 
criterion. We shall demonstrate that it is possible to overcome this limit by using 
entangled (correlated) photon number states, which have no classical analog. 

In classical interferometric lithography [2], when two coherent plane 
waves of laser radiation are made to intersect at an angle of 26, as shown in 
Fig. (1), interference fringes form with a spacing (pitch) of p = (A/2) sin 6. In the 
grazing incidence limit, 6 — > n/2, the minimum linear feature size that can be 

written is x mm = A/2. To see this, realize that the normalized exposure dose A at the 
substrate (proportional to the intensity) is given by the scaled interference pattern 
of two counter-propagating, grazing incidence, plane waves, A(x) = 1 + cos (2kx) = 

1 + cos 2(p. Here, k = 2tz/A is the optical wavenumber, <p = kx is the associated path- 
differential phase, and x the lateral dimension on the substrate. (We shall 
assume grazing incidence 6 — > n/2 for the rest of this work.) The Rayleigh 
criterion states that the minimal resolvable feature size occurs at a spacing 
corresponding to the distance between an intensity maximum and an adjacent 

intensity minimum [3]. The criterion then demands ^f™ = nil, from which we 

obtain x mm = A/2, as given above. This is the best resolution that can be achieved 
classically [2]. 

Recently, Yablonovitch and Vrijen (YV) have proposed utilizing "classical" 
two-photon exposure techniques to improve resolution [4]. The idea is that the 
uncorrelated (classical) two-photon absorption probability scales quadratically 
with the intensity [4]. Hence, the classical two-photon exposure dosage has the 

form, A c 2r = A 2 (x)/2= (l+cos2<p) 2 12 = f + cos2<p+ |cos4<p. This function has a 
term cos 4<p that oscillates in space with twice the frequency as the single-photon 

function, A{x). If the middle term containing the more slowly oscillating cos 2q> 
could somehow be eliminated, one would be left with a pattern resolution of 
x™ n = A / 4 ; a factor of two improvement. Since the number of elements writeable 

on a surface scales inverse-quadratically with the minimum feature dimension, 
this is an important advance. The approach of YV is to eliminate this middle 
term using a classical frequency modulation scheme. In the rest of this work, we 
show instead how to employ entangled photon number states to improve on this 
result. 
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It has been known for some time that entangled photon pairs, such as those 
generated by spontaneous parametric down conversion [4], have unusual 
resolving characteristics. This feature allows for sub-shot-noise interferometric 
phase measurement in theory [6, 7] and experiment [8]. In particular, Fonseca, et 
ah, recently observed the resolution of a two-slit diffraction pattern at half the 
Rayleigh limit in a coincidence-counting experiment [8]. What we will now show 
is that this type of effect is possible not only in coincidence experiments, but also in 
real two-photon imaging systems. Quantum entanglement is the resource that 
allows sub-diffraction-limited lithography. 

Consider the schematic set up for quantum interferometric lithography, 
illustrated in Fig. (1). We consider a two-port device with photons incident on a 
symmetric, lossless, beam splitter (BS) from the left in one or both ports A or B. 
The photons are then reflected by a mirror pair (M) onto the substrate (S) at the 
right. We can model the phase differential due to path-length differences between 
the upper and lower branches of the interferometer as a single phase shifter 

(PS) placed in the upper branch, which imparts a phase shift (p = kx. The two, 

photon paths converge on the imaging plane at an angle of off the horizontal 
axis, as show. We identify with the two input ports A and B, two photon annihila- 
tion operators d and b, respectively [7]. These operators obey the usual photon 
commutation rules, [d,d f ] = ^6,6 f =1 and d,b = 0. We can take the output 
electric field operator at the image plane on the right to be proportional to the sum 

of two output operators c and d from the upper and lower branches of the 
interferometer, C and D, respectively. Then the linear relationship between the 
two inputs and the two outputs can be expressed by a two-dimensional matrix 



equation, T 
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, where T is the input-output matrix. For our purposes, f 



can be thought of as being the product of matrices of the form, B = 
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, and P = 







which represent the unitary actions of the BS, M, 



and PS, respectively. Here, we have assumed a % phase shift on each reflection off 

a beam splitter or a mirror and a nil phase shift upon transmission through the 
BS, as required by reciprocity [9]. We assume that the phase differential between 
the two paths in Fig. (1) is represented in the single parameter <p of the PS. Hence, 
for Fig. (1), we have T = PRB. Relating input to output, we now deduce the output 
operators as c = (d - ib)e l(p I \[2 and d = (-ia + b)/ V2. Hence, the total, scaled, 
electric-field annihilation operator e, at S, is given by e = c + d= + e l9 ja + 

-j=(l-ie lf }b. Now, to compute the normalized, quantum, two-photon, exposure 

dosage A 2y , it is sufficient to compute the expectation values of moments of e and 
e\ since AT-photon absorption rates at the imaging surface will be proportional to 
the expectation values of S N = (e f ) (e)^ / AH the dosing operator. The quantum 
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theory of uncorrelated two-photon absorption was first worked out in 1931 by 
Maria Goppert-Mayer [10]. More recently, Javanainen and Gould [11] reported the 
equivalent theory for entangled two-photon absorption, while Perina, Saleh, and 

Teich [12] developed the entangled Af-photon theory. The expectation values of 8 N 
are taken with respect to input states that enter at the left of Fig. (1). The advan- 
tage of our operator-based approach is that all of the properties of the interferome- 
ter are encoded into the operator form of e = c + d, and one may then use different 
input states without having to recalculate the effects of the interferometer for each 
state [7]. 

Consider the input state = 1 1) A |0) B , where uncorrelated photons are 
incident one-at-a-time on the upper port A (Fig. 1). For this state we have a 
deposition rate of A lr (<p)= {Wi\^Wi) = l-sin2<p= 1+ cos(2<p + f ), which is the 
usual classical result, up to an unimportant phase shift. (Since the photons are 
uncorrelated, the interference pattern is the same as for |a) A |0) B , a classical 

coherent-state input [7].) Notice that A lr (<p) has a mean value of one, which 
normalizes the optical power deposited per unit length to unity. The classical two- 
photon deposition rate is then, A c 2y = (l+cos2<p) 2 12 = f + cos2<p+ \ cos 4<p, dropping 
the phase factor. We wish to use quantum interference to eliminate the slowly 
varying middle cos 2q> term. A simple choice of a nonclassical number-product 
state accomplishes this is \yf n ) = | 1) A | 1) B . This state is the natural output of a 
single-photon parametric down-conversion event [4]. The deposition rate for this 
state on a two-photon absorbing substrate has the form A q 2y = (i// n \8 2 \y/ n ) = 

l+cos4<p, as desired. The term containing the slowly oscillating cos 2<p has been 

eliminated, and we are left with only the cos 4<p term with the resolution of 

x^ m =A/4. In the inset of Fig. (1), we plot the classical one-photon pattern 

(dashed), the classical two-photon pattern (dotted), and our entangled two-photon 
(plain) results. Note that our quantum pattern has narrower features and half the 
peak-to-peak spacing of either classical curve. 

To understand the physics behind this improvement, as well as the role of 
entanglement, we adapt an argument used by Huelga, et al., in the context of 
atomic-clocks [13]. (A two-port Mach-Zehnder interferometer is isomorphic to a 
two-pulse Ramsey atomic-clock interferometer, under the SU(2) spinor rotation 
group [6].) Important to note is the form of the quantum state inside the interfer- 
ometer in Fig. (1), after BS but before PS, at the points A' and B'. Interference 
effects upon passage through a symmetric, lossless beamsplitter, cause the 
product number state \ys n ) = |1) A |1) B to become the entangled number state \y/ E ) = 

(|0) A ,|2) B , + |2) A ,|0) B ,)/ a/2. Here, the entanglement is between photon number and 

path; it is not possible to tell whether both photons took the lower path or both took 
the upper path [5]. This entangled state is sometimes called a "diphoton" or a 
"two-photon", and it behaves as single quantum-mechanical object, of photon- 
number two [14]. There are two indistinguishable paths this diphoton can take 
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through the interferometer, and so the amplitudes corresponding to these two 
paths will interfere. However, the diphoton will pass through the phase shifter 
only on the upper path, and this probability amplitude will acquire twice the 
phase shift as with a single photon process. Therefore, it is easy to see that the 

entangled state becomes, \yr E {(p))= (|0) c |2) fl + e 2i «'|2) c |0)J/ V2. This is the origin of 
the doubling of the resolution in the deposition rate A q 2y , given above. 

We can compute this same rate A| by a slightly more general argument. 
Suppose we ignore the process by which the two-photon entangled state is created 
in the interferometer, for the moment, and simply assume that the state \y/ E {(p)) 

has the maximally entangled form, given above, just to the right of PS. It now is 
not necessary to relate the output to the input operators. Instead, we can compute 

the deposition expectation function directly as (y^ E {(p)\S 2 \yf E {<p)) , using e = c + d, 

which gives the same result A q 2y , as before. It is now a simple matter to show that 

if an entangled number state of the form \y/ E {N)) = (\0) A ,\N) B , + \N) A ,\0) B .)/ -J2 is 

prepared inside the interferometer at the points A' and B' (Fig. 1), then the phase 
from the path differential accumulates N times as fast, producing \y/ E (N,(pf} = 

(\0) C \N) D + e m \N) c \0) D )l 42 at the output ports [13]. If the substrate is AT-photon 
absorbing, then the deposition function may be directly computed as, 

^M = {¥e{^M%We{^M) =l+cos2Aty . (1) 

All the slowly oscillating terms have been eliminated, leaving only the term 
cos 2N(p that has a resolution of AJ(2N). This is a remarkable factor of AT below the 
classical, Rayleigh limit of A/2. 

Of course, the natural question arises: how can one produce maximally 
entangled photon number states |y £ (Af)) of this form? For the case of AT =2, the 
answer is a parametric down-conversion event followed by a symmetric 
beamsplitter. For higher N, one approach is to use a material with a x m 
nonlinearity, in which a single photon is downconverted into N entangled 
daughter photons [12]. Another approach is to utilize a cascaded arrangement of 

N-l crystals with %< 2) nonlinearities [12, 15]. It is interesting to note that the 
natural output of a series of parametric downconversion events, or that of an 
optical parametric oscillator (OPO), has the form \N 1 2) A \N 1 2) B of a Hilbert space 
product [6,7]. In the case N = 2, only a simple linear beam splitter need be used to 
generated the maximally entangled form of \y/ E (N)y More generally, for N>2, 
one can construct such a state by operating on \N / 2) A \N / 2) B with a series otN, 
nonlinear, photon-photon, controlled-NOT gates [13, 16]. 

It is clear that the procedure described above is sufficient for etching a 
series of parallel lines with a factor of N reduction in the pitch. For more general- 
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purpose lithography, one would like to write arbitrary patterns in 2D. It is a 
simple matter to extend our technique to do this. Our approach makes use of the 
most general entangled state that can be generated by the above-mentioned 

techniques for fixed N, namely, \yr N>P {(p)}= (e lP,p \N - P) A \P) B + e i(N - p)<p \P) A \N - P)). 

With phase shifters, we can combine two such states interferometrically as 

|^iv,p,p'} = a N,p\VN,p) + Pn,p'\¥n,p') > with complex coefficients a and (3 such that 

2 i i2 

\a\ + \P\ = 1. This set of states can now be used as a basis to expand the target 

intensity pattern p(q>) in terms of the deposition components A^ p = 

(^n,p,p , \^n\ x ^n,p,p')- The task now is to calculate an appropriate set of these 

coefficients in order to fit a target function p((p). This can be done using a genetic 
algorithm optimization scheme [17]. For an example of this approach, let us 
approximate a ID square well "trench", given by the sample target function 

p((p) = 1- @{(p-7il2) + @(<p-37z:/2)where is the unit step function. The target 

function p((p) is depicted as a dark solid line in Fig. (2). The dotted line shows the 
best possible one-photon classical interferometric approximation. The solid curve 
is our quantum interferometric approximation, expanded with P e {1,...,5} while 
N = 10, and optimized with a genetic algorithm. Note that our entangled pattern 
(solid curve) is a much better approximation to the target trench function (dark 
solid line) than is the best classical approximation (dashed curve). It is a simple 
matter to extend this technique to 2D by adding orthogonal photon modes, as we 
shall show in a forthcoming work [17]. This method for producing arbitrary 
patterns in 2D is analogous to techniques used in commercial lithography — the 
Fourier components of the desired image are manipulated using phase masks, to 
insure that the actual image is as close to the target image as possible [18]. 

We make one final important note on the entangled AT-photon absorption 
process. Classically, the uncorrelated AT-photon absorption probability scales like 
I N , where / is the normalized classical intensity. This result is a consequence of 
the fact that the photons arrive independently of one another. Hence, the probabil- 
ity that the first photon arrives in an elemental absorption volume in space-time is 
proportional to /, and the probability a second photon will also happen to be in the 
same volume is also proportional to /, and so on, giving the I N law. For this 
reason, AT-photon absorption lithography with uncorrelated "classical" light is 
infeasible for high N, since extremely high flux densities are required [12]. This is 
not the case for AT-photon absorption with entangled photons. For two-photon 
absorption with entangled photon pairs, the absorption cross section scales as /, 
and not the I 2 that would be expected classically [11]. In general, this scaling law 
holds for arbitrary N. In other words, the AT-photon absorption cross-section, with 
N entangled photons, scales like / and not I N expected classically [12]. This result 
can be seen by the following heuristic argument: Recall that the photons are 
correlated in space and time, as well as number. Hence, if the optical system is 
aligned properly, the probability of the first photon arriving in a small absorptive 
volume of spacetime is proportional to /. However, the remaining N-l photons 
are constrained to arrive at the same place at the same time, and so each of their 
arrival probabilities is a constant, independent of I. Hence, although classical N- 
photon lithography requires unrealistically high optical powers, entangled N- 
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photon lithography requires only the same levels of power as the classical one- 
photon device. 

In conclusion, we have discussed the problem of entangled iV-photon 
absorption, as applied to interferometric optical lithography. We conclude that 
highly entangled photon states can be used to write features in an AT-photon 

absorbing resist which have a resolution of X/(2N). This result is a factor of N 
below the classical Rayleigh diffraction limit. Such states can easily be made for 
N=2 using optical parametric downconversion, and there are several possible 
approaches for implementing the scheme for higher N. We also indicate how the 
technique may by extended to write arbitrary 2D patterns. Entanglement turns out 
to be a useful resource, which can be employed in a technology such as lithogra- 
phy to overcome seemingly unbreakable constraints such as the diffraction limit. 
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FIGURE CAPTIONS 

Fig. 1 Interferometric lithography set up utilizing photons entering ports A 
and B. The photons strike the symmetric, lossless, beamsplitter (BS) and 
then reflect off the mirrors (M). The photon amplitude in the upper path 

accumulates a phase shift cp at the phase shifter (PS), before the two 
branches interfere on the substrate. Inset: The deposition rate A as a 
function of the phase shift q> for uncorrelated single-photon absorption 
(dashed), uncorrelated two-photon absorption (dotted), and entangled 
two-photon absorption (solid). Note that the classical two-photon curve 
has narrower features than the one-photon, but that the entangled two- 
photon has even narrower features still, and also shows the critical 
halving of the peak separation. 

Fig. 2 Here we show a target "trench" function (dark solid line), the best 
classical interferometric approximation (dashed), and our quantum- 
entangled interferometric approximation with ten entangled photons. 
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